Radiation 'damping' in atomic photonic crystals 
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The force exerted on a material by an incident beam of light is dependent upon the material's 
velocity in the laboratory frame of reference. This velocity dependence is known to be difficult to 
measure, as it is proportional to the incident optical power multiplied by the ratio of the material 
velocity to the speed of light. Here we show that this typically tiny effect is greatly amplified in 
multilayer systems composed of resonantly absorbing atoms (e.g. optically trapped ^^Rb), which 
may exhibit ultra-narrow photonic band gaps. The amplification of the effect is shown to be three 
orders of magnitude greater than previous estimates for conventional photonic- band-gap materials, 
and significant for material velocities of a few ms~^. 

PACS numbers: 42.50.Wk, 37.10.Vz, 42.70. Qs, 67.85.-d, 03.30.+p 



The force of radiation pressure is dependent upon the 
velocity of the body being pushed. For a perfectly reflect- 
ing mirror this velocity dependence appears as a kinetic 
friction term in the equations of motion for the mirror, 
and is due to the reduction in photon flux and frequency 
as observed in the material's rest frame, relative to the 
laboratory frame. This phenomenon was predicted some 
time ago by Braginski and Manukin [1 , where it was 
observed that the oscillatory motion of such a mirror 
connected to a wall via a spring would be damped in 
proportion to the power density of the incident beam. 
More recently, there was revived interest in this effect 
in reference to the precise interferometry experiments re- 
quired to detect gravitational waves (e.g. the LIGO and 
Virgo projects) [2 . However, as stated in [3 , to observe 
these velocity dependent terms, even in the case of a 
perfectly reflecting metallic mirror, the most favourable 
parameters lead to a laser power density so great that 
the mirror would be unlikely to remain intact. Therefore 
the question must be asked as to whether there are other 
physical systems where the fundamental velocity depen- 
dence of the force of radiation pressure could possibly be 
observed. 

Here we examine the radiation pressure experienced 
by a one-dimensional multilayered atomic structure [4 , 
with incident radiation of a frequency close to an atomic 
transition. This is done with the help of the Maxwell 
stress tensor [5 , which enables us to arrive at an exact 
expression for the pressure exerted by a light pulse, and in 
turn to asses the feasibility of measuring radiation damp- 
ing with such an atomic structure. For suitable choices 
of optical lattice period within this frequency window, 
an array of trapped ^^Rb atoms is known to exhibit an 
ultra-narrow photonic-band-gap [6 , with a width on the 
order of GHz [71 E] . Due to the very high sensitivity of 
the optical response (transparency or reflection) of such 



a multilayer to the frequency of the incident radiation, 
we observe that the velocity dependence of the force of 
radiation pressure is greatly enhanced, and can take ei- 
ther sign. This enhancement is three orders of magni- 
tude above that recently predicted for conventional one- 
dimensional photonic crystals [3 . We note that our re- 
sults are not specific to trapped ^^Rb, but apply equally 
well to any system exhibiting a photonic-band-gap on 
the same frequency scale. 

The rate of transfer of four-momentum, dP^^r^/dt^ to a 
medium which is possibly dispersive and absorbing, may 
be calculated from the electromagnetic fields in the vac- 
uum region outside of the medium, 



dt 
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where the surface element dSj points outward from the 
material surface [5 , and where the relevant components 
of the energy-momentum tensor are T^* = ceo{E x B)i 
andT^^' = eo[Sij{E^^c^B^)/2-EiEj-c^BiBj]. We start 
by considering a material planar slab at rest, upon which 
a beam of linearly polarized radiation of cross sectional 
area A is normally incident. Radiation of frequency uj en- 
ters the material through the surface at x = xi and exits 
through a similar surface at x = ^2 , with complex reflec- 
tion and transmission amplitudes r{uj) and t{uj). Aver- 
aging over a time interval At ^ uj~^ yields the average 
four-force experienced by the slab at rest, 



dt 



AeoEi ( [1 - ^(^) - n^)\ 
[l^R{uj)-T{uj)]xj 



(2) 



where R{uj) = |r(a;)p and T{uj) = |t(a;)p. Upon Lorentz 
transforming (pi), one obtains the corresponding expres- 
sion for the slab in motion with velocity V = Vx in the 



lab frame. In terms of the lab (primed) frame quantities, 

one has w = ^{1 -V/c)/{l + V/c)uj' = tjuj', Eq = t]E'q, 



and the average four-force becomes, 
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where P' = AcqcEq /2 is the incident radiation mean 
power. The force exact expression ([3| has a rather in- 
volved dependence on the velocity as it stems from R 
and T that depend on the velocity through the Doppler 
effect. A characteristic velocity range for cold atoms 
experiments is such that 77 :^ 1 — V/c^ so that the 
reflectivity and transmissivity can be expanded to the 
leading order in V/c around the lab frame frequency 
as, i?(V) - R{^') - {uj'V/c)R^^\uj'), and T(V) - 
T{uj') - {uj'V/c) T^^\uj'). This yields the simpler expres- 
sion on the right hand side of ([s]) where we set, 

R{uj') - T{uj') 

3R{uj') - T{u') - J [r^^^ {J) + T^i) (cjO) 

• R{yj') - T{uj') 

• 3R{u') - T{uj') + uj' (r^^^ {J) - T^i) (a;0) • 
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FIG. 1: Reflectivity and transmissivity for a multilay- 
ered atomic structures of 5.3797 x 10"^ unit cells with 
a=370.873xl0"^ m and 6=19.4807xl0"^m. Atomic pa- 
rameters are, ( ^'^Rb D2 line 5^ S1/2 -^ 5^ P3/2) 7e = 
27r X 6 X lO^Hz, cjq = c/Xo = 27r x 384.02 x 10^^ Hz and an 



For non-dispersive materials F^^^ and the the first three atomic density N/V = 6 x 10^^ m-^ ( A/" = 1.15 x 10 
terms in FA^. are the only contributions to the radiation 
pressure. These are positive and for a lossless medium of 
fixed reflectivity Rq reduce to the familiar results [TJ |2] 
respectively for the pressure force 2RqP' /c and the ra- 
diation pressure damping (friction) —ARqP'v/c^. The 
latter is typically much smaller than the former and 
for non-dispersive mirrors, where Rirjuo') ~ R{uj') and 
TiyTjuo') ~ T{uj')^ velocities as large as V ~ lO^ms"^ 
are needed to observe a few per cent velocity shift, cor- 
responding to a damping contribution to the force. For 
dispersive materials, on the other hand, we have an extra 
contribution to F\-. whose sign depends on the relative 
strength of the reflectivity and transmissivity gradients 
R^^\uj') and T^^\uj'). When R and T change signifi- 
cantly over a frequency range Acj <C uo' ^ this contribution 
could be substantial, and become the dominant term in 



F\-. whose sign, unlike other terms in the force expres- 
sion in (3), can then become negative or positive. This 
amounts, respectively, to either amplification or damping 
of the medium relevant dynamics [3[ [9] . Materials with 
optical reflectivity changes /S.R ~ 1 on the MHz range, 
i.e. (j^AR/Auj ~ 10^, would yield an appreciable velocity 
shift in the force even at V ~ ms~^. 

Periodic structures of trapped two-level atoms sepa- 
rated by vacuum fit the aforementioned parameter range 



quite well [7], and typical reflectivities and transmissiv- 
ities are plotted in Fig. [l] as a function of the scaled 
detuning S = {ujq — oj)/je from the atomic resonance 
transition cj^, where 7e denotes the excited state decay 
rate. These specific profiles correspond to a stack of al- 
ternating (complex) refractive indices [10] na{uj) c^ 1 and 
ni){uj) ~ Y^T^TSTrA/yp^^^Ty respectively with thicknesses 
a and 6, where JV = N/{VXl) denotes the scaled density 
of atoms. These multilayered atomic structures exhibit 
two pronounced photonic stop bands when the Bragg 
scattering frequency is not too far from the atomic tran- 
sition frequency. The first one develops from the polari- 
tonic stop band and has one edge at the atomic resonance 
frequency, the second one corresponds to the usual stop 
band of more familiar (non resonant) photonic crystals 
and develops around the Bragg frequency [7]. For the 
case of Fig.l the first and second stop bands lie respec- 
tively on the negative and positive detuning regions. At 
a given frequency u^ the radiation pressure is solely de- 
termined by the multi-structure optical response and Fig. 
^ shows both force components F.L and FA^. in the red- 
oetuned spectral region near the polaritonic stop-band in 
Fig. [1] In particular. Fig. ^p displays the force velocity- 



dependent contribution, FA^, which largely arises from 

the term proportional to R'^^^uj') —T'^^\uj'). This gradi- 
ent difference amounts precisely to 2xR^^\uj')-\-A^^\uj') 
and, unlike familiar dielectric photonic crystals struc- 
tures [11 , absorption dispersion A^^\uj') makes an ap- 
preciable contribution to the force at the edge(s) of the 
polaritonic stop-band where energy transfer to the sam- 
ple through absorption is largest [7 . Radiation pressure 
values can scale with V/c by as much as I^qJ = 1 x 10^ 
(ignoring the region very close to resonance) so that an 
appreciable ten per cent radiation pressure shift is pos- 
sible with velocities of ~ 3ms~^, that are now within 
experimental reach [12]. It is worth noting that the ve- 
locity contribution to the radiation pressure turns out to 
be at least three orders of magnitude larger than those 
obtained with traditional dielectric photonic crystals [3]. 
In addition, depending on whether the difference between 
R''^\uj') and T'^^^uj') is positive or negative alternating 
cooling and heating cycles become easily accessible in 
the appropriate spectral region [19^. We should note 




-2 xlO^ 



FIG. 2: Normalised x component of the force F as a 
function of detuning around resonance where radiation 
pressure force are strongest. The two panels show the 
force velocity independent {F}^^) and velocity dependent 

(^(\)) parts. 

that the above analysis is approximate due to the ex- 
pansion to first order in F/c, yet the above conclusions 
are not fundamentally altered if all orders of V/c are re- 
tained. This is shown in the following in the case of a light 
pulse, or sequence of pulses, of finite extent rather than a 
(monochromatic) plane light-wave. This is indeed to be 
examined to better asses the experimental feasibility and 



becomes furthermore important when opto-mechanical 
effects of radiation pressure associated with weak quan- 
tum light pulses are to be observed [H [131. The mo- 
mentum exchanged between a pulse of a specified extent 
and the multilayered atomic structure described above is 
here calculated by keeping all orders of V/c in the result- 
ing expression. At the position of the material surface, 
xi, the electric field of a pulse of central frequency Uc 
and half-width-half-maximum, C^/\n{i)^ can be writ- 
ten as a superposition of its incident and reflected parts 
as follows (the real part of this expression entering the 
energy-momentum tensor), 

E{t) = r ^i{u) (e^-(^i/^-^) +r(a;)e-^^(^^/^+^)) . 
Jo v27r ^ ^ 

(4) 
When the pulse is linearly polarized along the ?/-axis and 
the wavelength associated with the carrier frequency is 
much shorter than the extent of the pulse {uJcC./ c ^ 1), 
then a Gaussian pulse can be described as: ^(cj) c^ 

(/:(u;- 



/ioA^^^c>^/(^\/7r/2)ye-^^^^-^-)/2^) , where N is the 
average number of photons within the pulse. The trans- 
ferred momentum from a pulse incident on the structure 
entrance surface at Xi, as observed in the material's rest 
frame, is calculated as an integral over all time of the rel- 
evant integrated energy-momentum tensor component in 
(1). As an example, in T^^ one encounters the integral, 

/CO /"CO 

dt {E(tf + ^B{tf\ = / dw \^{u)f [1 + R{co)] 
-oo ^0 

(5) 
A similar contribution for the momentum transfer at the 
rear surface at X2, can be obtained starting from Q 
where we replace e^^(^i/^-^) +r(a;) e-^[^(^i/c+^)-^-(^) -^ 
t(c(;)e*'^^^2/c-t) rj^j^g ^Q^ integrated four-momentum im- 
parted by the pulse as seen in the slab rest frame can be 
calculated and in turn transformed as a four- vector into 
the lab frame. If the change in the medium velocity is 
negligible over the transit time of a single pulse, the net 
lab (primed) frame momentum transfer per pulse is given 

by, 



APf" 



eoAr] 



2^/1 - (V/c)^ Jo 



dLo'\$,{r]Lo') 



po' 



(6) 



where in the lab (primed) frame we have; C = C /rj] 
ujc = v^c^ V^' = I- R{W) - T{r]uj') + {y/c)(l + 
R{r]uj') - T{r]uj')); and V^' = 1 ^ R{r]uo') - T{r]uo') + 
{y/c){l — Rirjuo') — Tirjuo')). This recovers well know 
results for a transparent lossless slab [T4| in the limit 
for which V/c -^ and 77 ^ 1. The plane wave limit 
in ([3| also emerges namely when £ ^ 00, if AP^' is 
first divided by the proper time interval for a square 
pulse containing the same amount of momentum, and of 
the same maximum amplitude as the incident Gaussian 
pulse, Ar = -s/tt /2{C/c). The results presented in figure 
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FIG. 3: Normalised impulse, AP' • x acquired from a 
3m long Gaussian pulse impinging upon a multilayered 
atomic structure. The spectral region has been divided in 
two parts (a) and (b) for clarity. The detuning Sc between 
the pulse carrier frequency (ujc) and atomic resonance 
(cjo) are in units of 7e. All other parameters are as in fig. 

in 







[s] illustrate the momentum transfer, AP • x^ in units of 
the total momentum contained in the incident pulse in 
the lab frame, Nfik'^. They substantiate the findings of 
the previous section: Fig. [3^ shows a difference in the ra- 
diation pressure force per pulse ~ 10% {5c ^18) relative 
to a stationary structure (V = 0) while a similar differ- 
ence of ~ 15% is shown to occur away from resonance, 
where (5c ^ 65. Radiation damping effects from multilay- 
ered atomic structures forming Bragg mirrors are best in- 
vestigated in optical lattices created by the interference 
of multiple laser beams, which cool and localize atoms 
at the lattice sites [15]. For a sufficiently long trap the 
(Bragg) mirror may be envisaged as made of an array of 
disks spaced by half the wavelength A of the confining op- 
tical lattice U{x) and filled with atoms in the vibrational 
ground state of the lattice wells. The disks thickness d is 
essentially given by the rms position spread around the 
minima of the potential U{x) while the transverse size 
D ^ ^J^No(l|l\Lpd is determined by the in-well atomic 
density p and the number of atoms Nq loaded into a trap 
of length L [20 . For typical densities p ~ VS^^cvrT^ 
and filling factors hja c^ 0.1 a transverse size D in the 
range 35 ^ 110 /im is obtained for Nq ~ 10^ ^ 10^ when 
L ^ \ cm. This sets the incident beam waist Wq and in 
turn its Rayleigh range xr = Trw^/Xo that should clearly 
compare with the overall mirror length L. Taking waists 
Wo ^ 20 ^ 90 /im yields xr ^ 0.25 -^ 3.2 cm, which can 



safely be made larger than L for the cloud with more 
atoms. Ultracold samples of Na atoms, e.g.^ may fur- 
ther suit the case since No ^ 10^ associated with an even 
shorter transition wavelength A^ can be attained [16 . 

After loading the atoms into a ID optical lattice, whose 
counterpropagating beams are taken to be far (red)- 
detuned from atomic resonance ^Tj , the lattice is further 
set into motion dragging along the atoms. Such a mo- 
tion may be achieved by changing the relative frequency 
detuning, Auj of the two laser beams, which corresponds 
to a lattice velocity v = Auj /2k where k is the aver- 
age wavenumber. Through such a scheme velocities on 
the ms~^ range can be reached [12]. Upon passage of a 
light pulse, the transferred momentum AP' induces co- 
herent oscillations in the center of mass of the atomic 
wavepacket in the lattice wells. This in turn induces a 
periodic redistribution of the power difference between 
the two counterpropagating lattice beams that can be 
measured [15]. In particular, every half a cycle of the 
atomic oscillations, the momentum of the atomic ensem- 
ble changes by 2 AP' and this corresponds to a coherent 
scattering of a number of photons ANph from one of the 
two counterpropagating laser beams of the optical lat- 
tice to the other which is given by hk ANph = |AP'|. 
Of course, by sending a train of Nip light pulses equally 
spaced by the period of the atomic harmonic oscillations 
so to always push the atoms at the right time, ANph 
would be amplified by the factor Nip. Ail-optically tai- 
lorable light pressure damping and amplification effects 
can here be attained in the absence of a cavity, which 
makes multilayered atomic structures not only interest- 
ing in their own right but also amenable to a new opto- 
mechanical regime. The large per-photon pressures that 
can be observed compare in fact with state of the art 
micro |T7] and nano ^18] optomechanical resonators, yet 
involving (atomic) masses that are various orders of mag- 
nitude smaller. 
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